The Race Between Stars and Quasars in Reionizing Cosmic Hydrogen 
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The cosmological background of ionizing radiation has been dominated by quasars once the Uni- 
verse aged by ~ 2 billion years. At earlier times (redshifts z > 3), the observed abundance of bright 
quasars declines sharply, implying that cosmic hydrogen was reionized by stars instead. Here, we 
explain the physical origin of the transition between the dominance of stars and quasars as a generic 
feature of structure formation in the concordance ACDM cosmology. At early times, the fraction 
of baryons in galaxies grows faster than the maximum (Eddington-limited) growth rate possible for 
quasars. As a result, quasars were not able to catch up with the rapid early growth of stellar mass 
in their host galaxies. 



Introduction. Observations of the Lycv forest and the 
quasar luminosity function at redshifts z < 6, suggest 
that quasars dominated the production rate of hydrogen- 
ionizing photons only after the Universe has aged by 1-2 
billion years [see Fig. 5 in Ref. [l( . The sharp decline in 
the observed comoving density of bright quasars at red- 
shifts z > 3 implies that supermassive black holes (BHs) 
could not have produced sufficient UV photons to reion- 
ize cosmic hydrogen by a redshift of z re ion = 10.9 ± 1.4, 
as required by the WMAP5 data on the microwave back- 
ground anisotropies 0]. It is therefore widely believed 
that stars have reionized cosmic hydrogen [l|, Here 
we suggest a simple physical explanation for this phe- 
nomenological inference. 

Preliminaries. Quasars are powered by the accretion 
of gas onto massive BHs 0, @| ■ The remnant BHs they 
leave behind are observed in the nuclei of all bulged galax- 
ies at the present time @ . BH growth is accelerated dur- 
ing episodes of galaxy mergers, when cold gas is infused 
to galactic nuclei by tides; this association is predicted by 
computer simulations @, B @] and supported by observa- 
tions • The observed correlation between the mass of 
nuclear BHs and the depth of the gravitational potential 
well of their host galactic bulges (as inferred from the ve- 
locity dispersion of stars) [llL su ggests that their growth 

Merely ~ 5% of the 



was ultimately self-regulated [12|, 
energy output from a bright quasar is sufficient to unbind 
the cold gas reservoir from its host galaxy d, [l]| • 

Bright quasars are inferred to be radiating close to 
their limiting Eddington [lij ] luminosity, Le- Estimates 
of their BH masses imply typical luminosities L in the 
range of - (0.1-1) x L E [see Fig. 6 in Ref. 15]. The 
Eddington limit is derived by equating the repulsive ra- 
diation force on an ionized gas element to the attractive 
gravitational force on it towards the BH [l||, 
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where <jt is the Thomson cross-section for electron scat- 
tering and M is the BH mass. In a spherical geometry, 
the radiation force associated with L > Le would gener- 
ate an outflow and inhibit accretion onto the BH. Given 



an unlimited fuel reservoir, the growth rate of a BH would 
be regulated by the maximum luminosity that allows it 
to accrete gas, L E . 

The luminosity of a quasar is related to the gas ac- 
cretion rate M a cc through L = eAf acc c 2 , where e is the 
efficiency for converting the rest mass of the accreting 
gas to radiation. At high accretion rates, the BH lumi- 
nosity is expected to approach its limiting value Le [13] , 
as inferred for bright quasars [H|]. Theoretical models 
imply that at L > 0.5Le, the accretion disk is puffed-up 
by radiation pressure 3, 1|| and its geometry enters the 
quasi-spherical regime for which the Eddington limit is 
derived. 

The growth rate of the BH mass, M = (1 — e)M acc is 
related to the Salpeter [J] time, 
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For constant values of e and L/Le, the growth of the 
black hole from a seed mass M see d over time At is expo- 
nential, with M = Mgeed exp(At/tfi). Interestingly, the 
growth time tE has no explicit dependence on BH mass, 
and so it remains the same irrespective of whether mul- 
tiple seeds grow in parallel or they coalesce to grow as a 
single BH. 

At the high accretion rates of bright quasars, the cool- 
ing time of the gas is shorter than its accretion time 17] 
and a thin accretion disk with a high radiative efficiency 
[20I ] forms. The existence of such disks has been con- 
firmed recently by microlcnsing observations [21.] ■ The 
radiative efficiency of the gas is dictated by the inner 
boundary of the disk at the innermost stable circular or- 
bit [2^1, from where the gas plunges into the BH and 
its viscous dissipation rate diminishes For a non- 

spinning BH, e = 5.7%, while for a maximally-spinning 
BH, e = 42% 0, During the prodigious growth 

phase of quasars, the BH is expected to be spun-up 
quickly by the infalling gas, but the orientation of the 
angular momentum vector of the disk might vary consid- 
erably between different accretion episodes. Therefore, 
we adopt an intermediate value between these efficien- 
cies in our fiducial numerical example. The fact that 
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the cosmic BH mass budget grew through a luminous 
accretion mode of a high radiative efficiency (e ~ 10%) 
and not through a hidden mode of a low radiative effi- 
ciency, is demonstrated by comparing the radiation en- 
ergy density produced by quasars to the BH mass density 
in the local Universe [23J, |24| . A high radiative efficiency 
is also implied by the clustering and abundance statistics 
of quasars j25|. 

Stars form through fragmentation of cold gas in disks 
of galaxies [26] . The growth of the stellar mass budget oc- 
curs on the dynamical timescale of the host galactic disk 
[42I ]. The most vigorous mode of star formation (star- 
burst activity) is also realized in gas-rich galaxy mergers 
within which cold gas is concentrated into a compact re- 
gion jjjj]. The fragmentation (Jeans) mass is lowered to 
the scale of stars only in environments that are denser by 
> 4 orders of magnitude than the mean cosmic density p; 
hence, the associated dynamical time is guaranteed to be 
shorter by > 2 orders of magnitude than the age of the 
Universe [~ (Gp)" 1 / 2 ] at all redshifts. The undelayed 
growth in the stellar content of galaxies through merger- 
driven starbursts is ultimately limited by the global as- 
sembly rate of cold gas into galaxies. This rate is in turn 
dictated by the cosmological build-up of dark matter ha- 
los in which gas may cool. 

The deposition of energy and momentum by stars or 
quasars into their gaseous environment may expel gas 
from the host galaxy and further regulate the growth of 
the stellar and BH mass budgets. However, here we focus 
on the maximum allowed growth in both components and 
show that the mass density of stars could have grown 
faster than that of quasar BHs at very high redshifts. 
There are indications that these maximum rates were 
realized in nature. The existence of bright quasars with 
BH masses M > 1O 9 M at z ~ 6 [27], when the age 
of the Universe was only tnub = 10 9 yr[(l + z)/7]~ 3 / 2 , 
implies that the exponential (Eddington-limited) growth 
of their BHs must have persisted over a time interval At 
that covers much of their cosmic history Q , 
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Despite this persistent growth (43[, the observed quasar 
population did not supply sufficient UV photons to reion- 
ize the Universe at z > 6 even though it dominates the 
production of ionizing photons at z < 3 [l[ . This implies 
that the stellar mass budget was able to grow faster than 
quasars at early cosmic times relative to their respec- 
tive growth at late times. Below we show that this phe- 
nomenon is a generic outcome of the evolution of struc- 
ture in the standard ACDM cosmology. 

The production rate of ionizing photons is proportional 
to the growth rate in the mass densities of BHs and stars 
per comoving volume in the Universe [3]. While the 
minimum growth time of BH mass is i^, the shortest 



growth time of the stellar mass budget, t ga \, is dictated 
by the assembly rate of cold gas into galaxies. At late 
cosmic times, tE is much shorter than t gal , leading to 
a feedback regulated mode of BH growth in which the 
supply of cold gas is a limiting factor 12, HI • Below 
we show that at early cosmic times the situation was 
reversed. 

Mass assembly rate of galaxies. The minimum mass 
of galaxies in which BHs and stars form is dictated by 
cooling considerations. At the redshifts of interest here, 
the cooling time of the gas is shorter than its dynamical 
time for dark matter halos that have a virial temperature 
T V ii > 10 4 K, above the threshold for atomic hydrogen 
cooling [2q . Galaxies above this threshold are believed 
to have hosted the bulk of the sources that have reionized 
the Universe. Early on, the minimum virial temperature 
of galaxies might have been reduced by an order of mag- 
nitude through molecular hydrogen cooling, but H2 can 
be easily dissociated by UV photons • After reioniza- 
tion, the minimum T v ; r is expected to have increased to 
~ 10 5 K, owing to photo-ionization heating of the inter- 
galactic medium from where galaxies are assembled [3(| ■ 
The fraction of baryons available to make stars equals the 
mass fraction of dark matter that virialized in halos above 
the minimum T v j r . We denote this fraction by f(z) and 
calculate the timescale for the assembly of cold gas into 
galaxies from the time derivative / = (df / dz) / (dt / dz) , 
where dt/dz = 8.4 x 10 7 yr[(l + z)/10]~ 5 / 2 , so that 
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In the Press-Schechter formalism [31], / = 
erfc[(5 c (z)/ v / 2CT(M min )], where S c = 1.686/£>(z) is 
the collapse threshold for an overdensity linearly- 
extrapolated to z = 0, D(z) is the linear growth factor of 
density perturbations (with D(0) — 1), and cx(M m i n ) is 
the root-mean-square amplitude of linearly-extrapolated 
density perturbations at z = on the scale of a sphere 
from which the minimum galaxy mass is assembled. 
Refinements to the halo mass function 32j have a 
negligible effect (< 15%) on our results for t ga \. 

Results. Figure 1 shows the ratio (t ga \/tE) as a func- 
tion of redshift z for the concordance ACDM cosmology 
2]. The galaxy growth time is calculated as the frac- 
tion of matter that is incorporated in dark matter ha- 
los with virial temperatures above 10 3 K (dashed line), 
10 4 K (solid), and 10 5 K (dotted). We adopt [e/(l - e)] = 
32%(L/Le), corresponding to a BH spin-averaged radia- 
tive efficiency of e = (6 + 42)/2 = 24% for L = L E or to 
e = 5.7% for the characteristic value 15j of L/Le ~ 0.2. 
This choice is conservatively smaller by a factor of ~ 2 
than the lower limit implied by observational data on 
high- redshift quasars [251 ]: a corresponding increase in 
this value by some factor would have lowered the plotted 
curves by the same factor and strengthened our conclu- 
sions. The three vertical lines mark the central value 
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Redshift 

FIG. 1: Redshift evolution of the growth time (in tE units) 
for the fraction of matter that is incorporated in dark mat- 
ter halos with virial temperatures above 10 3 K (dashed line), 
10 4 K (solid), and 10 5 K (dotted). Conservatively, we adopt a 
value of [e/(l — e)] = 0.32(L/L_b), which is a factor of ~ 2 
smaller than the lower limit implied by observations of high 
redshift quasars [25l ] . The vertical lines correspond to the cen- 
tral value (solid line) plus or minus one standard deviation 
(dashed lines) for the redshift of reionization, z rc i on , based 
on the WMAP5 data [2j. Prior to reionization, the assem- 
bly of cold gas into galaxies occurs at a faster rate than the 
Eddington-limited growth of BHs. 



(solid line) plus or minus one standard deviation (68% 
confidence; dashed lines) for the redshift of instantaneous 
reionization, as inferred from the WMAP5 data 0]. Prior 
to reionization, the assembly of cold gas in galaxies oc- 
curs at a faster rate than the Eddington-limited growth 
of BHs for all three threshold values of T vlr . 

At early cosmic times, galaxies are on the exponen- 
tial tail of the Press-Schechter mass function, and the 
mass of cold gas in them grows faster than tE- BH 
growth through gas accretion is not able to track this 
early rapid growth. At late times, quasar growth is lim- 
ited by the much slower rate at which fresh cold gas 
is infused into the centers of galaxies through episodic 
mergers. For the post-reionization case of T VK > 10 5 K, 
we find that < 0.5i ga i only at z < 6. It takes sev- 
eral e-folding times for quasars to build-up their mass 
density and catch up with the UV production by star 
formation in common galaxies at lower redshifts. De- 
spite the initial delay in their growth, quasars eventu- 
ally dominate the production rate of ionizing photons 
by z ~ 3 1]. Within 1-2 billion years after the big 
bang, the growth in the cosmic mass density of BHs 
starts to be regulated by feedback on their gaseous en- 



vironment (owing to their increased mass) and by the 
depletion of their cold gas reservoir (aided by the declin- 
ing merger rate of galaxies). The characteristic time for 
doubling the mass of galaxies traces the age of the Uni- 
verse (~ 2.3 x 10 9 yr[(l + z)/4]~ 3 / 2 ), and keeps increas- 
ing relative to 4e with decrea sing redshift. The mode of 
feedback-regulated growth [H, ll2l.[l3T] becomes critical at 
these low redshifts. 

During the early history, when BH growth is not yet 
regulated by feedback or by the exhaustion of the cold gas 
reservoir in galactic nuclei, the comoving mass density of 
accreting BHs (pbh) grows at a rate, 

PBH = PBli/tE + Ascod, (5) 

where p SO cd is the formation rate density of BH seeds. 
If the BH seeds cannot ionize the Universe on their own 
(i.e., if reionization did not result from the accretion lumi- 
nosity of stellar-mass BHs) , then Fig. 1 shows that stars 
are required and able (ZBj to make reionization happen 
by z = 10.9 ± 1.4. 

At late times, the BH mass density grows exponentially 
until feedback and the consumption of cold gas (which 
were not included in Eq. O start to regulate its growth 
[46| . By requiring that quasars match the UV production 
rate per comoving volume of stars only as late as z ~ 6-3 
Pl, we infer that they fall short of matching it at z ~ 11 
by orders of magnitude. The transition redshift between 
the early domination by stars and the late domination 
by quasars is sensitive to uncertain parameters (which 
may also be redshift dependent), such as the unknown 
feedback strength and the radiative efficiency of these 
source populations. 

The radiative efficiency of stars and BHs is obviously 
different. The total number of emitted ionizing photons 
per baryon incorporated into stars ranges [36| between 
4 x 10 3 [for a present-day (Pop II) mass function] and 
~ 10 5 [for metal- free massive stars (Pop III)]. Supermas- 
sive BHs produce ~ 2x 10 7 e ionizing photons per accreted 
baryon [371 ]. The observation that quasars dominate the 
cosmic UV production rate only at low redshifts [1[ im- 
plies that their BH formation efficiency out of cold gas 
in galaxies is lower than that of stars by several orders of 
magnitude. This is likely caused by angular momentum 
which distributes the cold galactic gas in a large-scale 
disk that fragments into stars long before tidal or viscous 
torques enable a small fraction of this gas to feed the cen- 
tral BH. The inference of a low BH formation efficiency 
is confirmed by data on the cumulative mass bud gets of 
supermassive BHs and stars in the local Universe [35| . 

The fundamental limitation presented by Fig. 1 ap- 
plies also to an early population of stellar-mass black 
holes, because those would still require extensive accre- 
tion in order to reionize the Universe. The comoving 
luminosity density of any BH population is simply pro- 
portional to its cumulative accretion rate, and the Ed- 
dington limitation on the accretion time is independent 
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of the mass distribution of the early BHs. Indeed, the 
smallness of the unresolved component of the soft X-ray 
background places severe constraints on the cumulative 
mass density of an early population of accreting BHs ir- 
respective of its mass distribution [4(| . 

Observable signatures. The delay associated with the 
Eddington- limited growth of black holes in the early Uni- 
verse can be probed through a number of observational 
methods. First, gravitational wave signals from coalesc- 
ing BH binaries at high redshifts - which are detectable 
by LISA [13 and possibly also by Advanced LIGO [Hj], 
can be used to probe the mass function of massive BHs as 
a function of redshift. Second, future X-ray missions such 
as the International X-ray Observatory [48[ , and infrared 
telescopes such as JWST (49[, or new ground-based in- 
struments, may detect fainter quasars at higher red- 
shifts than those accessible with present-day telescopes. 
Third, BHs that escaped from their dwarf galaxies by 
gravitational-wave recoil at early cosmic times, are po- 
tentially detectable in the Milky- Way halo through the 
compact star clusters that surround them (39[. Finally, 
extending current simulations of quasar growth 0, 0] to 
the higher redshifts of reionization (z > 12) would refine 
predictions for the ionized bubble si zes, which may be 
probed by future 21cm observatories [371 kl| . 
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